Absence of skew scattering in two-dimensional systems: 
Testing the origins of the anomalous Hall effect 
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We study the anomalous Hall conductivity in spin-polarized, asymmetrically confined two- 
dimensional electron and hole systems, focusing on skew-scattering contributions to the transport. 
We find that the skew scattering, principally responsible for the extrinsic contribution to the anoma- 
lous Hall effect, vanishes for the two-dimensional electron system if both chiral Rashba subbands 
are partially occupied, and vanishes always for the two-dimensional hole gas studied here, regardless 
of the band filling. Our prediction can be tested with the proposed coplanar two-dimensional elec- 
tron/hole gas device and can be used as a benchmark to understand the crossover from the intrisic 
to the extrinsic anomalous Hall effect. 
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Introduction. — The observed Hall resistance of a mag- 
netic film contains the ordinary Hall response to the ex- 
ternal magnetic field and the anomalous Hall response 
to the internal magnetization. Although the anomalous 
Hall effect (AHE) has been used for decades as a ba- 
sic characterization tool for ferromagnets, its origin is 
still being debated, also in the context of a closely re- 
lated novel phenomenon, the spin Hall effect [l|, 0, 01 ■ 
Three mechanisms giving rise to AHE conductivity have 
been identified: (1) an intrinsic mechanism based solely 
on the topological properties of the Bloch states orig- 
inating from the spin-orbit-coupled electronic structure 
, (2) a skew-scattering mechanism originating from the 
asymmetry of the scattering rate and (3) a side-jump 
contribution, which semiclassically is viewed as a side- 
step-type of scattering and contributes to a net current 
perpendicular to the initial momentum Q. 

Recent experimental and theoretical studies of tran- 
sition-metal ferromagnets and of less conventional sys- 
tems, such as diluted magnetic semiconductors, oxide 
and spinel ferromagnets, etc., have collected numerous 
examples of the intrinsic AHE and of the transition to the 
extrinsic AHE dominated by disorder scattering [s^ . The 
unambiguous determination of the origin of the AHE in 
these experimental systems is hindered, in part, by their 
complex band structures, which has motivated studies of 
simpler model Hamiltonians, such as the two-dimensional 
(2D) Rashba and Dirac band models 0, [l^, El H [ll ■ 
Attempts to describe all the contributions to the AHE 
within the same framework have yielded farraginous re- 
sults, however. So far a rigorous connection of the more 
intuitive semiclassical transport treatment with the more 



systematic diagramatic treatment, providing a clear-cut 
interpretation of the intrinsic, skew, and side-jump AHE 
terms, has only been demonstrated for the Dirac Hamil- 



tonian model |14| . 



In this Letter we calculate the transport coefficients 
in these two complementary approaches for asymmetri- 
cally confined 2D electron and hole gases in the pres- 
ence of spin-independent disorder, finding perfect agree- 
ment. The motivation for the study of these systems 
is threefold: First, they can be represented by simple 
spin-orbit-coupled bands, which, similar to the Dirac 
Hamiltonian model, allows us to unambiguously identify 
the individual AHE contributions. Second, the extrinsic 
skew-scattering term vanishes for a two-subband occu- 
pation in the case of the Rashba 2D electron gas and 
for any band occupation for the studied 2D hole gas. 
This provides a clean test of the intrinsic AHE mech- 
anism and of the transition between the intrinsic and 
skew-scattering-dominated AHE. Finally, we propose a 
2D electron gas/2D hole gas coplanar magneto-optical 
device in which the unique AHE phenomenology found 
in our theoretical models can be systematically explored 
experimentally. 

Model Hamiltonians. — We study the following 2D 
model Hamiltonians: 
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with m being the effective in-plane mass, <Tj the 2x2 
Pauli matrices, k± — kx i: iky, h the exchange field, 
an the spin-orbit coupling parameter, and V{r) a spin- 
independent disorder potential. The exponent rt = 1 (3) 
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describes a 2D electron (hole) gas [Tsl] . The eigenenergies 
of the clean system are E± = h^k"^ /2m± yJK^ + (Q!„fc")2. 
The eigenvectors in the clean system take the form 

k (cos 0, sin <j)) and 



= exp(ik • r)|uj) with k 



1 / ±ie-™"*\/A±7^ 



(2) 



where A — ^Jh? + [unk'^Y . We now define k±{E) as the 
wave number for the ± band at a given energy E and de- 
fine X± = X{k±). If E is not specified, it is assumed to be 
the Fermi energy. We consider the model of randomly lo- 
cated (5-function scatterers, V{r) = J^i ^(i" ~ R-i) with 
Rj random and disorder averages satisfying {Vi)dia = 0, 
{V?)di. = + 0, and (F/)dis = + 0. This model 
is different from the standard white-noise disorder with 
(l^k'kP)dis = Tit^Q ^ where is the impurity concentra- 
tion and other correlators are either zero or related to this 
correlator by Wick's theorem. The deviation from white 
noise in our model is quantified by V\ ^ 0, and is neces- 
sary to capture part of the skew-scattering contribution 
to the AHE. 

Semiclassical approach. — We sketch here the semiclas- 
sical procedure used in the calculation, for further details 
we refer to Ref. O. The multi-band Boltzmann equation 
in a weak electric field E is given by 

dh^^^__df^ (3) 



dt 



± is the subband index, and 



J:^, J d^k'/{27:f Loiv ifi - fv) is the impu- 



where I 
/[/Icon = 

rity collision integral. The distribution function // is 
the sum of the equilibrium function and a correction, 
/; = fi + gi . The scattering rates ujw are related to the 
T-matrix elements through loui = 27r/?i |T;/;p5(e;/ — e;), 
where Tin = {^'\y\^i), and are eigenstates of the 
complete Hamiltonian, and |/) = of the disorder- 

free Hamiltonian. 

Skew scattering. — Skew scattering appears in the 
Boltzmann equation through the asymmetric part of the 
scattering rate, i.e., low ^ ujin fil. The scattering rates to 
second and third order in disorder strength are given by 
ujiv = 4'^ +4^^+- • • , where4'^ = 2^?* |')dis(5(e,- 
£;') is symmetric. Here Vin — {V\V\l). We break up the 
third-order contribution into symmetric and antisymmet- 
ric parts. We ignore the first, since only the second gives 
rise to skew scattering. This antisymmetric term is given 
by [iil 



,(3a) 
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(4) 

The solution of the Boltzmann equation ([3]) is found by 
first looking at the deviation of the distribution function 
from equilibrium (l^ . 
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Assuming that the transverse conductivity is much 
smaller than the longitudinal one (A^ ^ B^) and sub- 
stituting Eq. ([5]) into Eq. ^ one finds A^^ = rj and 
= where 
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For symmetric Fermi surfaces, the skew-scattering con- 
tribution to the conductivity tensor at zero temperature 
can now be expressed using the scattering times. 
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The calculation of (r^) ^ and (tJ 



(8) 

uses the matrix 



elements of Eq. Q . To simplify the notation we define 
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Jo 
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where all momenta are taken on the Fermi surface. Note 
that in Eq. ^ the magnitude of k" can be different from 
that of k' or k since the Fermi momenta of different bands 
do not coincide. 

The matrix elements appearing in Eq. ([9]) can be cal- 
culated directly from the basis functions, yielding 
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from which we obtain 
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where is related to the density of states of each band 
at the Fermi energy, iy^)^^ — T? jm ± n(a„fcJ~"^)^/A±. 
The symmetric part of the scattering rates to second 

(2) 

order in the disorder potential is given by ujly ~ 

2t: /hniV^\{u'^\u'^,)\^5{ei — ev). The relaxation times are 
found by inserting this into Eq. ^ and Eq. into 
Eq. ([7]). For n — 1, i.e., for the 2D electron gas, the 
relaxation rates are then 
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FIG. 1: (color online). Diagramatic representation of the 
skew-scattering contribution to Uyx- Both current vertices, 
denoted by squares, are renormalized by ladder vertex cor- 
rections. 



where p, = —fi. If both subbands are occupied, the last 
factor in Eq. (fT3)) vanishes and there is no skew-scattering 
contribution. If only the majority subband is occupied 
{Ep < h), {t^)~^ is non-zero and skew scattering con- 
tributes. For the skew-scattering Hall conductivity and 
the longitudinal conductivity we obtain in this case 
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If n = 3, i.e., for the 2D hole gas, we obtain 
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and skew scattering vanishes irrespective of band filling. 

Microscopic approach. — Within the diagramatic Kubo 
formalism the skew-scattering contribution to the off- 
diagonal conductivity is obtained from the expression 
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where the bare velocity vertex factors in the linear-in-k 
Rashba model are given by 
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As shown in a previous study ^14] , the skew-scattering 
contribution proportional to Vi/{niV^) corresponds to 
the diagrams shown in Fig. [1] where the current ver- 
tices jx , jy on both sides are the bare velocities Vx , Vy 
renormalized by ladder vertex corrections. Only the 
skew-scattering diagrams with a single third-order ver- 
tex, shown in Fig.[l] contribute to order Vi/{niVQ). All 
other terms from a ladder-type summation of third-order 
vertices are smaller because they are either not of the or- 
der l/rii or of higher order in Vi/Vq. The sum of the 
skew-scattering vertices (i.e., the bold/red part of Fig. [T]) 



gives 
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In the linear Rashba model we find i/+/A+ ~ v-fX-, 
implying that skew scattering vanishes if both subbands 
are occupied. In the case that only one subband is occu- 
pied the evaluation of Fig. [T]to order Vi/{niVQ) yields 
exactly the same expression for a^"" as in the semiclas- 
sical Eq. (fTSjl . The only effect of the ladder vertex cor- 
rections is to rcnormalize each bare velocity by a factor 
of 2(^,2 + A2_)/(3/i2 + Xl) which reduces to a factor of 1 
in the limit of small aikp and to a factor of 2 in the limit 
of small h. 

For the 2D hole-gas model Hamiltonian ([1]) with n ~ 3 
the bare velocity vertex factors are 
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Here the vertex corrections disappear because integrals 
of the type J2k ^k'^x,yG^ = vanish. This impHes the 
absence of skew scattering for any subband filling [l^, 
consistent with the semiclassical result. We note that the 
same consistency between semiclassical and microscopic 
quantum theory calculations for the studied 2D models is 
also obtained for the intrinsic and side-jump terms simi- 
lar to the results in the graphene model 14|; the longer 
details of those calculations will be shown elsewhere and 
are in general agreement with Ref. [l^. 

The abscence of the skew scattering is akin but not 
equivalent to the results of spin-Hall-effect calculations 
in 2D systems 17[. For the Rashba 2D electron gas the 
disappearance of the DC spin Hall conductivity is guar- 
anteed by sum rules that relate the spin current to the 
dynamics of the induced spin polarization [3, [l^ . In the 
case of a charge current no similar sum rule is known. As 
we have shown, the skew-scattering contribution in fact 
becomes finite when the minority band is depleated. The 
vanishing of the Hall conductivity in the Rashba 2D elec- 
tron gas for Ep > h is attributed to the simplicity of the 
Hamiltonian. In particular the relation Z/+/A+ = i^_/A_ 
does not hold generally beyond the case of the linear-in- 
k Rashba coupling. The abscence of skew scattering in 
the 2D hole system has a different origin: Due to the 
cubic dependence of spin-orbit coupling on momentum, 
the matrix elements, Eq. ()10|) . in the antisymmetric part 
of the collision term behave like sin(30 — Scj)')- Together 
with the sm{(f) — cj)') dependence of the velocity factor in 
Eq. ([7]), this makes the integral over k' vanish. 

Our results predict that the AHE in 2D electron and 
hole systems can be dominated by contributions indepen- 
dent of the impurity concentration, for which the anoma- 
lous Hall resistance is oc a"^. We also predict that in 
the Rashba 2D electron gas with only one subband oc- 
cupied the extrinsic skew-scattering contribution, lead- 
ing to anomalous Hall resistance proportional to cr"^^ , is 
non-zero. Note that this term has not been identified in 
previous works that considered only white-noise disorder 
M Hi Since its corresponding conductiv- 



ity contribution is inversely proportional to the impurity 
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and spin-orbit coupling in order to control the transi- 
tion between the intrinsic and extrinsic AHE regimes. 
The exploration of the AHE in the 2D hole gas is more 
challenging due to the expected sub-micron spin diffu- 
sion length in this system but may still be feasible in the 
proposed experimental setup. 
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FIG. 2: (color online). Top panel: Top- view schematics of 
the Hall bar with coplanar 2D hole and electron gases. Spin- 
polarized carriers are generated by shining circularly polar- 
ized light on the p-n junction. Center bottom panel: Cross 
section of the heterostructure containing p-type and n-type 
AlGaAs/GaAs single junctions. The left band diagram cor- 
responds to the unetched part of the wafer with the 2D hole 
gas, the right band diagram shows the 2D electron gas in the 
etched section of the wafer. 

concentration, the skew-scattering mechanism can dom- 
inate in clean samples. 

Proposed experimental setup. — The unique phenome- 
nology of the AHE in the studied 2D systems, in partic- 
ular the sudden disappearance of skew scattering when 
the Fermi level crosses the depletion point of the minority 
2D Rashba band, represents an opportunity for a clean 
test of the presence of intrinsic and extrinsic sources of 
the AHE and of the transition between these two regimes. 
In the absence of 2D ferromagnetic system with Rashba 
like spin-orbit interatciton, we proposed an experimental 
setup for this test as shown in Fig. [S] The device is based 
on a AlGaAs/GaAs heterostructure containing a copla- 
nar 2D hole gas/2D electron gas p-n junction. The cross 
section of the heterostructure and corresponding band 
diagrams are shown in the lower panels of Fig. [2] (for 
more details see Ref. H). Under a forward bias the junc- 
tion was successfully utilized as a light-emitting-diode 
spin detector for the spin Hall effect [3|. Here we pro- 
pose to operate the junction in the reverse-bias mode, 
while shining monochromatic, circularly polarized light 
of tuneable wavelength on the p-n junction. The pho- 
togenerated spin-polarized holes and electrons will prop- 
agate in opposite directions through the respective 2D 
hole and electron channels. The longitudinal voltage and 
the generated anomalous Hall voltage can be detected by 
the successive sets of Hall probes, as shown in the upper 
panel of Fig. O For the 2D electron gas the macroscopic 
spin diffusion length allows to use standard lithography 
for defining the Hall probes. Surface or back gates in 
close proximity to the 2D electron system can be used 
to modify the effective 2D confinements, carrier density. 
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